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Abstract. We consider the trigonometric Felderhof model, of free fermions 
in an external field, on a finite lattice with domain wall boundary conditions. 
. . . The vertex weights are functions of rapidities and external fields. 

' We obtain a determinant expression for the partition function in the special 

' case where the dependence on the rapidities is eliminated, but for general 

^ ~ ^ external field variables. This determinant can be evaluated in product form. 

In the homogeneous limit, it is proportional to a 2-Toda r function. 
» I ' Next, we use the algebraic Bethe ansatz factorized basis to obtain a product 

' expression for the partition function in the general case with dependence on 

all variables. 



^ ' 0. Introduction 

04 I In [1], Korcpin introduced domain wall boundary conditions for the six vertex 

' model on a finite square lattice, and proposed recursion relations that determine 

I the corresponding domain wall partition function. In [2], Izcrgin obtained a de- 

' terminant solution of Korepin's recursion relations. At the free fermion point, the 

, six vertex domain wall partition function can be evaluated explicitly in product 

' form [3]. In the homogeneous limit, it is proportional to a 1-Toda t function [4]. 

(~| I In this work, we look for analogous results in the context of the trigonometric 

CIh' limit of Felderhof 's model [5], which is a vertex model of free fermions in an external 

field. In section 1, we recall the definition of the model in the parametrization of 
Deguchi and Akutsu [6], and formulate it on an x lattice. There are four sets 

eof complex variables: horizontal and vertical rapidities {ui, Vj}, and horizontal and 
vertical external field variables {a^, where G {1, 2, ... , The weight 

^ I Wij of the vertex Vij at the intersection of the i-th horizontal line and j-th vertical 

line depends on the difference of the rapidities, Ui — Vj, but depends on the external 
^ , fields, ai and I3j, separately. 

' In section 2, we impose domain wall boundary conditions and obtain an Izergin- 

type determinant expression for the domain wall partition function, under the re- 
striction that the difference of any two rapidity variables is a multiple of 27r^— 1, 
but for general {ai,/3j}. This expression can be evaluated in product form. In the 
homogeneous limit, it is proportional to a 2-Toda r function [7,8]. 

In section 3, we use the factorized basis of the algebraic Bethe ansatz, [9, 10], 
to obtain a product expression for the domain wall partition function for general 
{ui,Vj} and {ai,(3j}. 

0.1. Abbreviations. In the rest of this paper, DWBC stands for 'domain wall 
boundary conditions' and DWPF stands for 'domain wall partition function' . Zj^^ 



X 



■'TF 

I N X N lattice. 

the DWPF of the six vertex model on an x iV lattice. 



is the DWPF of the trigonometric Felderhof model on an x lattice. Zg^^ is 
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^TFres -^TF ^ With restrictions on the rapidities as in Equation[3] -^^Fres hom 
is the homogeneous version of Z^p^^^ with all horizontal external field variables 
equal, and all vertical external field variables equal. 



1. The TRIGONOMETRIC FELDERHOF MODEL 

1.1. The lattice. We work on a square lattice consisting of N horizontal and N 
vertical lines. We label the horizontal lines from top to bottom, and the vertical 
lines from left to right. 

We assign the i-th horizontal line an orientation from left to right, a complex 
rapidity variable Ui and a complex external field variable a^. We assign the j-th 
vertical line an orientation from bottom to top, a complex rapidity variable Vj and 
a complex external field variable (3j. 



Ml, ai 

""2, "2 



MJV, CtN 



Vi V2 

f3i 132 



VN 

I3n 



Figure 1. An N x N square lattice, with oriented lines. To each line, we attach 
two complex variables, a rapidity and an external field. 



1.2. Vertices. Each line intersects with N other lines. A line segment between 
two intersections is a bond. To each bond, we assign a state variable, namely an 
arrow that points along the orientation of that line segment or against it. The 
intersection of the i-th horizontal line and the j-th vertical line, together with the 
four bonds adjacent to it, and the set of arrows on these bonds, is a vertex Vij. 




^ ^ ^ 

ai 62 C2 

Figure 2. The non-zero weight vertices of the trigonometric Felderhof 
model. The black arrows indicate state variables. The white indicate 
line orientation. Notice which vertex is ai and which is 02. 
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1.3. Weights. To each vertex Vij, we assign a weight Wij, that depends on 1. 
The orientations of the four arrows on the bonds of that vertex, 2. The difference 
of rapidity variables flowing through the vertex, and 3. The two external field 
variables flowing through the vertex. To satisfy the Yang-Baxter equations, only 
the six vertices shown in Figure [2] have non-zero weights [6]. In the parametrization 
of [6] , the non-zero weights are 



(1) ai{ai,(3j,Ui,Vj) = 1- aiPjc"' bi{ai, f3j,Ui,Vj) = ai - (3je'^' 

a2iai,l3j,Ui,Vj) = e"'""^ b2{ai, l3j,Ui,Vj) = l3j - aiC^'''"^ , 

In the sequel, we will drop the dependence on the variables, when that is clear 
from the indices. Unlike the six vertex model, the vertex weights of the trigono- 
metric Felderhof model are not invariant under reversing the directions of all the 
arrows. 

1.4. DWBC. As in the six vertex model, the DWBC are such that all arrows on 
the left and right boundaries point inwards, and all arrows on the upper and lower 
boundaries point outwards. 



Figure 3. Domain wall boundary conditions 

1.5. DWPF. Given a 2-state vertex model, such as the six vertex model or the 
trigonometric Felderhof model, the DWPF on an x lattice, -^j^^^c defined 
as the sum over all weighted configurations that satisfy DWBC. The weight of each 
configuration is the product of the weights of the vertices 



(2) 



'DWBC 



E n ' 

config- ^vertices 



It is also possible to define DWBC and DWPF's in vertex models with more 
state variables [11-13]. 



2. A DETERMINANT FORM OF THE RESTRICTED DWPF 

2.1. The Korepin-Izergin procedure. As a first step towards computing Z^p^ , 
we follow the Korepin-Izergin procedure: 1. Specify a set of properties that fully 
determine Z^P^ , 2. Conjecture a determinant expression for the required Z^p^ , 
and 3. Show that the conjectured expression satisfies the required properties. 

It turns out that it is not obvious how to follow the above procedure for general 
values of all variables. The reason is that an Izergin-type determinant solution is 
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tightly related to the Korepin-type properties, one of which is that the DWPF is 
symmetric under permuting the variables on any two parallel lattice lines. 

In the six vertex model, and in models discussed in [11-13], this condition is 
automatically satisfied because the vertex weights are invariant under reversing the 
directions of all arrows, so the two a-type vertices, which are involved in proving 
this symmetry, have the same weight. In the trigonometric Felderhof model, there 
is no such invariance for general values of all variables, and we need to impose 
restrictions on at least some of the variables. 

Our plan is to restrict the variables to a point where the Korepin-Izergin pre- 
scription works. We claim that there is no Izergin-typc determinant expression for 
Z^P^ , for general values of rapidities and external fields. 



2.2. Restrictions. We require 



(3) 

The restrictions in Equation [3] are satisfied by choosing the difference between 
any two rapidity variables to be a multiple of 27r\/— 1, or equivalently, by simply 
setting all rapidities to zero. The external field variables remain free. By eliminating 
the dependence on the rapidities, the weights are now much simpler and can be 
written as 



(4) ao,ij = ai(Q!i,/3j, 0, 0) = a2(ai, /3j, 0, 0) = 1 - 0;^/?^, 

= 6i(ai,/?j,0, 0) = -62(ai,/3j, 0,0) = Q!i - /3j, 

co,^J = ci(a„/3j,0,0) = C2(a„/3j,0,0) = ^1 -af^l -/^^ 
Under these conditions, z!^p^ becomes Z^p^^s- 



2.3. Korepin-type properties. Z^p^^^ is fully determined by the following prop- 
erties 

1. It is symmetric in the elements of each of the sets {a}, and {/?}. 



2. It is a polynomial of degree {N — 1) in a;, up to a factor of ^Jl — of, where 
1 < i < N , and in f3j, up to a factor of ^1 — /?|, where 1 < j < A^. 

3. It satisfies the recursion relation 



(5) 



^NxN 
^TF,ret 



N 




f \ 

N 








1=1 




Vii^n ) 



^(W-l)x(Ar-l) 

''TF,res,{'mn) 



AN-l)x{N-l) 



for any m,n e {1, . . • , N}. The subscripts (mn) in -^r^ _^es,(mn) 
the variables that are not present in the reduced partition function 
4. It satisfies the initial condition Zp^\^^ = cq.h. 



indicate 



2.4. Izergin-type determinant solution. Using the notation ajy] — ai — aj, 
etc, the properties in subsection 12.31 are satisfied by 
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(6) Z, 



NxN 
TF,res 



det 



Yll<i<j<N 



L (10,ijbo.ij ) I 



<iJ<N 



<i<j<N 



{a^ - aj){Pj - Pi) 



n 

I, l<k<N 



X det 



(it 



aiPj){ai - Pj) 



1 i<i,i 



<N 



2.5. Remarks on proof of Equation[6l The proof proceeds along the same hnes 
as Izergin's proof, which is discussed in detail in the literature, including [12, 17] 
and references therein, it suffices to outline it. Since the four Korepin-properties in 
[m fully determine ^^f^,, all we need to do is to show that the right hand side 
of Equation [6] satisfies each of these properties. 

Properties 1 and 2 can be checked precisely the same way as in the case of the 
six vertex model [12, 17]. Property 4 can be checked by inspection. Property 3 can 
be checked as follows. 

Expanding the determinant in the right hand side of Equation 6 along the first 
row, we obtain 



N 



(7) Z, 



NxN 
TF,re; 



Co,l) 



n^2 "[ij] llj<i Pm UJ>^ Pm 



^iN-l)x{N-l) 



Because of the DWBC, the vertex at the upper right corner must be either bi or 
ci. By choosing ai =Pn, we eliminate the possibility of a hi vertex, and restrict the 
allowed configurations as follows 1. The upper right corner is a ci(q!i, /Jat, 0, 0) = 
COjiAT vertex, 2. The right most column, apart from the upper right vertex, is a 
set of ai{ai, Pn,0,0) = ao,iN vertices, where {2 < i < N}, and 3. The top row, 
apart from the upper right vertex, is a set of a2{ai, Pj,0, 0) — ao,ij vertices, where 
{l<J<iV-l}. 

Letting ai = Pn in Equation [7l we obtain 



(8) 



^NxN 
^TF,ret 



C0,1N —Tr „^f-1 ^ 



nN YlN-1 Q 

j=2"[lj] lli=l P[Nj\ 



y{N-l)x(N~l) 
^TF,res,(lN) 





'N-1 




' N 


Co,lN 


n "O'lj 




n "-o-jN 











7(JV-l)x(Ar-l) 
^TF,res,{lN) 



which is the recursion relation of Equation [S] as we expect. Thus the determinant 
expression on the right hand side of Equation [6] satisfies all Korepin-properties. 

2.6. Further check on Equation [6l The proof of Equation [6] outlined above 
made use of a recursion relation obtained by freezing the upper right corner. We 
could have also chosen to freeze the upper left corner. Let us check that Equation 
[6] satisfies that second recurrence relation as well. 

Because of the DWBC, the vertex at the upper left corner must be either a2 or 
ci. By choosing aiPi =1, we eliminate the possibility of an 02 vertex, and restrict 
the allowed configurations as follows 1. The upper left corner is a ci{ai, Pi,0, 0) = 
Co, 11 vertex, 2. The left most column, apart from the upper left vertex is a set 
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of 62(0^, 0, 0) = — &o,ii vertices, where {2 < i < TV}, and 3. The top row, 
apart from the upper left vertex, is a set of bi{ai, f3j, 0,0) = &o,ij vertices, where 
{2 < J < N}. 

Letting ai/3i = 1 in EquationO we obtain 



(9) 



^NxN 
^TF,ret 



ai/3i = l 



C0,11 



(11^12 ao,ij^o,ii] (11^12 aojii^oji ] 

n^2 "[ij] n^2 Ail] 



Co, 11 



3- 
■ N 

n ^o.ij 



' N 



-,(Ar-l)x(Ar-l) 
^TF,res, (11) 

^{N-l)x{N-l) 
^TF.res, (11) 



which is what we expect. Equation[9]is a second recursion relation for Z^p^^^, and 
provides an independent check of Equation [6l 



2.7. On a determinant form with more general parameters. It is natural 
to look for a determinant expression for the DWPF with less restrictions on the 
rapidities than in Equation [Sj We were unable to find any such expression, even 
for the simplest variations on the conditions of Equation [3l such as allowing only 
one rapidity, such as ui, to be free, and so forth. This, of course, is not a proof 
that no such generalization exists, but only that, if there is one, it is unlikely to be 
of the Izergin form of Equation (6] 



2.8. The homogeneous limit. In the homogeneous limit ^ a, and /3j — > /3, a 
standard procedure gives 



(10) Z,^^^ 



TF,res,ho7n 

N(N-l) 



(-1)- 



((a-/3)(l-a/3)) [Vl-a^V^^] 



N 



det 



[[da) [dp) [ [a - P)il - a(3) ) ^ 



l<i,j<N 



2.9. 2-Toda r-function. Because the determinant in Equation llOl is bi-Wronskian, 
with partial derivatives in two complex variables, it is straightforward to show [8], 
using the Jacobi identity for determinants, that it is a r-function of the 2-Toda 
partial differential equation 



d tn+iTn-1 

As mentioned earlier, the homogeneous limit of Izergin's determinant expression 
of Zq^^ is proportional to a bi-Wronskian with partial derivatives in one complex 
variable, and therefore is a r-function of the 1-Toda partial differential equation [14]. 

This observation was used in [14] to study the free energy of the six vertex model 
in the presence of DWBC. Since Z^p^ can be computed explicitly in product 
form using the algebraic Bethe ansatz, as we will see below, the free energy can 
also be computed explicitly, and the relationship with 2-Toda remains a curious 
observation. 
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2.10. On enumeration. As is well known, Z^-^ can be used to enumerate al- 
ternating sign matrices (ASM's) [15, 16]. At the free fermion point, Zq-,J^^ 2- 
enumerates ASM's [16]. 

Z^p^ can also be used to enumerate ASM's, but because the model is yet again 
a free fermion model, one can easily show that here too one obtains 2-enumerations. 

2.11. A product form for Z^p^^^. The determinant in Equation [6] factorizes 



(12) det (X^^^) 



<i<j<N 



(1 - ai(3j){ai - (3j) 



The simplest way to see this is to notice that there is a change of variables that 
allows one to re- write the determinant in Equation [6] in Cauchy forn:0. 

This is reminiscent of the factorization of Izergin's determinant in the six vertex 
model, at the free fermion point [3]. We attribute the factorization of Equation [12] 
to the fact that the trigonometric Fclderhof model is a free fermion model. 

From Equation [6] and Equation 1121 we obtain 



(13) 



yNy.N 

-'TF,res 



n ^ 

l<i,j<N 



The simple form of the factorized result in Equation 1131 suggests that a similar 
result may hold in the general case with dependence on all parameters. This will 
be the topic of the next section. 



3. Product form for general Z, 



NxN 
TF 



Unlike Z^F 



7NxN ryNxN invarfant under permuting adjacent variables. All 

expressions in this section are valid only for the ordering shown in Figure [l] Ex- 
pressions of Z^p^ with different orderings are related by factors of vertex weights. 

3.1. Definitions. Consider weights which depend only on the vertical variables 
{(3} and {v} 



a2,^J = e^^-^^ - ft/3, 



(14) 



C2,ij 



1-/3^/1-/3 e 



The i?-matrix for the trigonometric Felderhof model is 



n h. 

Rij{Pi,Pj,Vi,Vj) 



ai,y U U 

ci^ij 

C2,ij b2,ij 

0,2, ij , 
The monodromy matrix, T0.1...7V, for N sites, is 



(15) 



7o,i...Af(a,M) = RoNia, (3n,u,vn) ■ ■ ■ Rniia, Pi,u,vi) 

^ ( Ai...jv(q!,m) Bi...jv(a,u) ^ 
[ Ci...jv(a,w) Di,,,N{a,u) )^ 



^This remark is due to Ch Krattenthaler. 



8 A CARADOC, O FODA, M WHEELER AND M ZUPARIC 

As explained in [17], Z^p^ can be expressed in terms of the creation operators, 
Bi...N, as 

(16) Z^^"" = (l|Bi...Ar(ai,ui)...Bi...^(a^,7/w)|0) 

where 

(ii = ®u ( 1 ) , = ®f=i (^0 ' = ( ) . ^ "^^-i (^0 

The expression in Equation [16] is not easy to evaluate directly, since the creation 
operators are sums containing 2^ terms, and each term is a tensor product acting 
in all of the spaces 1, . . . , A^. 

3.2. A factorizing matrix. Following [9] , we define an initial factorizing F-matrix, 
Fi^2...N, by 

Fi^2...n{Pi; (32, ■ ■ ■ , 13n]Vi;v2, . . ■ ,vn) = e^"^ + ep^^Ti,2...Ar(/3i, i^i) 

^^^^ " ( G2...n{Pi,Vi) D2...n{Pi,Vi) ) ^ 



From that, the full i^-niatrix, i^i...Ar, is defined recursively by 



(18) Fi...Ar(/3i, . . . ,/3iv;wi, . . . ,UAr) = F2...nFi^2...n = 

■ ■ ■ = F(jv-i)Ar-F'(Ar-2),(Ar-i)Ar ■ • • Fi^2...N 

where 



(19) F(Ar_i)Ar(^(Ar_i),^7v;W(Ar-l),f7v) = 

r 1 







1 



C2,(N-1)N b2,(N~l)N 



I 















a2,(N~l)N 



{N-l)N 



3.3. A twisted monodromy matrix. The full i^-matrix, Fi. jv, is now used to 
construct a twisted monodromy matrix 



(20) ToA...N 



Fi...nTo^i,,,n ^a; (3i, . . . , Pn;u;vi, . . . ,vn^ ^iT."^jv 



3.4. Twisted creation operators. We define twisted versions of the Bethe ansatz 
operators as follows. ^i...Ar(a,u) — Fi..jvv4i..jv(a, V' Using the nota- 
tion ai.oj = 1 — ct(3j e"""^, 02, oj = e"""^ — a(3j, etc, where the label indicates 
dependence on the horizontal variables {a, u}, one can show that 



(21) 



N , 

(8) 

,7=1 ^ 









Q2,jQQ2,Qj 



DWBC TRIGONOMETRIC FELDERHOF MODEL 



9 



(22) Bi,„N{a,u) 
^ ( h 



E(8) 

1=1 j<i 



( 0] 

[ cm J 



3>l 



Q Q2,ZjQ2,0jQ2,jf 



(23) Ci...N{a,u)-- 

E(8) 

1=1 j<i 



02,1 





' 


c2,o; 




u 





J 


(8) 











62.13 

Q a2,0j 

*^2,ij 



(24) i?i...Ar(a,M) 



3.5. Remarks on proof of Equations I21|, \22\ 1231 and 1241 One first verifies by 
direct computation that the above formulas hold for ^12(0;, u), Bi2{a, u), Ci2(a, u) 
and Di2(a,M). This becomes the basis for a proof by induction, in which each 
formula is proven individually. For example, to prove the formula for the twisted 
operator i3i...7v, we observe that, by construction, the untwisted operator -Bi...Ar 
satisfies 



(25) Bi...jv(a,w) = ^2...Ar(a, w) 



r ] r 62.01 ] 

n + B2...N{a,u) ' 

Multiplying Equation 1251 from the left by i^i...Ar, and from the right by F^^j^, 
we find that i3i...jv(a, m) satisfies 

(26) Bi,„Nia,u)Fi,„NF-\ = 

Fi...^F-:^I^A2.Ma,u)\^ ^^^^ ]^ + B2...A.(a,^) ( Y a2,oi]} 

From Equation 1181 we have Fi^^^^F^^ff = ^2...Ar^i,2...Ar^2~ V' Equation 
becomes the following matrix equation in space 1 

(27) Bi,,,N{a,u) I ~ ] ~ ^ I = 

L C2...n[Pi,vi) D2...n[Pi,vi) )^ 

[ C2...Nif3i,vi) D2...N{f3i,vi) ] ^ 

{^--("'")(cl )/^--("'") ( Y all]. 

This is a recursion relation for i3i.. jv(a, u) in terms of twisted operators over 
the TV — 1 spaces 2, . . . ,N. We use it to prove Equation [22] inductively, as follows. 

First, we postulate that the expressions of Equations \21\ \22\ 1231 and |24I hold 
over the — 1 spaces 2, . . . ,N. Next, we substitute them into Equation 1271 From 
that, the expression in Equation 1221 for i?i...7v is seen to be the unique solution to 
Equation 1271 Repeating this procedure for the other twisted operators, we prove 
the postulate over the N spaces 1, . . . , A^. 
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3.6. A recursion relation. Since (l|Fi...jv — IlKfc '^2jfc(l| and Vl'^) ~ l'^)' 

■Wx, 
TF 



lj<k 

we can re- write Z^p^ in terms of the twisted creation operators 



(28) 



zNxN (l|-Bl...Ar(Q;i, Ui) . . . Bi,,,n (aN , Un)\0) 



Following [10], we use the expression in Equation 1281 together with the explicit 
expression for the twisted i3-operator, to derive the recursion relation 



AT N-l , 

(29) z^-p^ = c,,m n n n ^^.^o n «M^4F:f^-r^'^ 

1=1 j"=l k=^i 2,ifc h>i 

where the subscripts in (^at^)^ indicate the omission of the variables {ajv, mjv} 
and {Pi, Vi}. 

3.7. Remarks on proof of Equation [29l Acting with _Bi.. jv(aAf , "A') on |0), in 
Equation 1281 we immediately find 



N 



(30) z!^^'' = — 



n "2 



^ [ij] ^i...Af(ai,Wi)...Bi...7v(a(w-i),M(Ar-i)) [i*] [tj] 

j=i 

Using the identities ^J,i j ct;" ^ j,; j = (ct^)^ = Oj one obtains 

(31) [ii] i?i...Ar(ai,Ui)...i3i...Jv(a(Ar-l),ti(Ar-l)) [ii] = 



'w-1 










n "2j* 




■J-j- a2,iA; 














1 fc<i 



-Bi...(i_i)(i+i)...Ar(ai, Ui) . . . i?i...(i-i)(i+i)...Ar(a(7V-l)7 ""(A-l)) 



Equation 1291 comes from substituting Equation 1311 in Equation 1301 and using 
the fact that 



(Af-l)x(Af-l) 
TF, (Ni) 



1) I\j<iO'2,jil\i^kO''^, 



(32) Zrpp - rr ^ 

^ ' [lj<k<^2,jk 

[-ij] ^i...(i-i)(i+i)...Af(cH''"i) • ■ • -^i...(»-i)(»+i)...Jv(a(Ar-i),'"(7v_i)) l^tj ] 



3.8. Solution of recursion equation. It can be shown that the following product 
expression satisfies Equation 1291 
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(33) Z^^^ 

N 



ne'=(---Vl-«|Jl-/3| 



k fc=i 



yi<j<k<N 



N 



I. k=l 



l<j<k<N 



3.9. Remarks on proof of Equation 1331 From Equation 1331 we have 



(64,) Z^rj^p^ (JY^J 



N 

n ' 

k=i+l 



N-1 



1 



Using Equation 1341 and after considerable manipulation, one recovers 



r^NxN 



N N-1 ^ 

k^i 



AT 

E 



Af-1 



n 



a2{aj,Pi,Uj,Vi) 

L a2(aj,aAr,Mj,UAr) 



n 



h {Pk,Pi,Vk,Vi) 



yNxN 



Finally, we observe that 



(36) 



N 

E 



N-l 

a2{aj,Pi,Uj,Vi) 



Yl_b2{Pk,aN,Vk,UN) 



Jl b2{f3j,l3k,Vj,Vk) 



Jl b2{f3-i,(3k,v.i,Vk) 



'7V-1 

which can be checked by noticing that both sides are polynomials of degree N — 1 
in the variable a^, and furthermore the equality is satisfied at the N points a^v = 
/Jje"^""™, where j — 1, . . . ,N. Equation 1361 means that 



N 



( N-1 



(37) E 



T-r a2{aj, f3i,Uj,Vi) 
-IJ- a2(Q;j,aAr,Uj,UAr) 



n 



b2iPk,aN,Vk,UN) 
b2{Pk,l3i,Vk,Vi) 



= 1 



which, when substituted in Equation I35[ causes it to collapse to the recursion 
relation Equation 1291 as required. 
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